Diverging equilibration times in long-range quantum spin models 
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The approach to equilibrium is studied for long-range quantum Ising models where the interaction 
strength decays like r~ a at large distances r with an exponent a not exceeding the lattice dimension. 
For a large class of observables and initial states, the time evolution of expectation values can be 
calculated. We prove analytically that, at a given instant of time t and for sufficiently large system 
size N, the expectation value of some observable {A)(t) will practically be unchanged from its initial 
value (A)(0). This finding implies that, for large enough N, equilibration effectively occurs on a 
time scale beyond the experimentally accessible one and will not be observed in practice. 



Equilibration is one of the central concepts of ther- 
modynamics, but our understanding of the underlying 
microscopic processes is still far from complete. Stud- 
ies of the approach to equilibrium can be traced back to 
Boltzmann's work [l[ in the early days of statistical me- 
chanics, and they are also closely related to the micro- 
scopic foundations of the second law of thermodynamics. 
Modern facets of this topic, attracting much attention 
in the recent research literature, are relaxation after a 
uantum quench Q , thermalization close to integrability 
u , or typicality as a foundation of quantum statistical 
mechanics to name but a few. 

Although experience shows that equilibration takes 
place in the vast majority of situations, exceptions have 
always attracted particular interest. A famous example 
is the numerical investigation of the time evolution of a 
chain of nonlinearly coupled oscillators by Fermi, Pasta, 
and Ulam [f| where recurrent behavior was observed, but 
no sign of equilibration. Another particularly interesting 
case, going under the name of quasistationary states, has 
triggered intense activity, reviewed in [f|, in the field of 
long-range interacting systems. Long-range refers here 
to interactions decaying at large distances rasa power 
law r~ a with a positive exponent a not exceeding the 
spatial dimension d of the system. Gravitating masses or 
Coulomb forces in the absence of screening are important 
examples of long-range interactions. The term quasista- 
tionary is used to describe metastable states whose life- 
time diverges with increasing system size N. The physi- 
cal importance of quasistationary behavior should be ob- 
vious: For a sufficiently large system, the transition from 
a quasistationary state to equilibrium takes place on a 
time scale that is larger than the experimentally accessi- 
ble observation time. Hence, equilibrium properties will 
not be observed, and instead the statistical properties of 
the quasistationary regime are of interest. 

Most of the studies of quasistationary states have fo- 
cused on the Hamiltonian Mean-Field model Q, a toy 
model consisting of classical XY spins (or plane rota- 
tors), each coupled to every other at equal strength. The 
typical scenario observed is that, for a suitable class of 
initial conditions, the total magnetization of this spin 



model rapidly relaxes to some quasistationary value dif- 
ferent from its equilibrium value, remains there for a very 
long time, and finally decays to equilibrium (see Fig. 33 
in [6]). More recently, quasistationary states have been 
observed in gravitational sheet models Q and their ex- 
istence has been argued to be generic for a large class 
of classical long-range systems [9[. Virtually all finite- 
system results in the field have been obtained by numer- 
ical techniques, supplemented by analytical calculations 
in the N — > oo Vlasov continuum limit [f|. 

In this Letter, analytic results on equilibration in long- 
range quantum spin models are reported. Surprisingly, 
the class of models studied permits exact analytic results 
for both, finite and infinite systems, even beyond what 
has been achieved in the classical case. In principle, the 
time evolution is found to be recurrent for any finite sys- 
tem size N, and no equilibration occurs in a strict sense. 
But since the recurrence times diverge exponentially with 
iV", they will be irrelevant for large enough systems and 
effective equilibration is observed. Similar to the qua- 
sistationarity in classical systems, we find that the sys- 
tem remains close to its initial state for extremely long 
times and effective equilibration times diverge with in- 
creasing system size. This is, to the best of the author's 
knowledge, the first observation of this kind of behav- 
ior in quantum systems, and it prepares the ground for 
further studies, including the above-mentioned quantum 
quenches or inquiries into the foundations of statistical 
physics of long-range systems. The results complement 
the peculiarities of long-range quantum spin systems in 
equilibrium reported in 10]. 

Groundbreaking advances have been witnessed in re- 
cent years regarding the experimental realization of spin 
models, where ultracold gases in optical lattices have 
been used to engineer spin systems whose parameters 
can be tuned. Long-range interactions between the 
constituents can be created by either trapping dipolar 
molecules 11| , or by inducing a nonpermanent dipole mo- 
ment when shining appropriately tuned laser light onto 
atoms [l2| . Future experiments of this kind might permit 
the experimental observation of diverging equilibration 
times in long-range quantum spin systems. 
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The model. — The main result of this Letter can be 
proved for a large class of (generalized) Ising models with 
long-range pair interactions, defined on various types of 
lattices and for arbitrary lattice dimensions. To keep the 
presentation simple, we discuss only the one-dimensional 
case explicitly. Higher-dimensional lattices and other 
generalizations are discussed below. 

Consider N identical spin-1/2 particles, attached to 
the sites {1,...,N} of a finite one-dimensional lattice 
with periodic boundary conditions. The corresponding 
quantum dynamics takes place on the Hilbert space 

N 

w = ®c?, (1) 

i=l 

where the <Cf are identical replicas of the two-dimensional 
Hilbert space of a single spin-1/2 particle. The time evo- 
lution on H is generated by the Hamiltonian 

jV N/2 N 

H N =M N Y,J2 e C7X*f+j - h (2) 

i— 1 j — 1 i—1 

with a z denoting the z component of the Pauli spin op- 
erator and ft, G R an external magnetic field in the z di- 
rection. The index i + j is to be considered modulo N to 
account for the periodic boundary conditions. The e(j) 
are pair coupling constants depending on the distance j 
of two spin operators on the lattice, and we assume 

lim e(j) = 0. (3) 

A typical example we have in mind is the Dyson model 
13] with algebraically decaying couplings, e(j) = j~ a , 
but, in contrast to Dyson's work, exponents a smaller 
than 1 are also allowed. More general classes of interac- 
tions will be discussed below. For exponents a > 1, the 
interaction is absolutely summable, 

oo 

X>c?)i<°°, ( 4 ) 

3=1 

and this case has already been treated in EMU. Here 
we are interested in what is sometimes called strong long- 
range interactions with < a < 1 where the sum in Q 
diverges. In this case, the double sum in © gives rise 
to an infinite energy per spin {Hn)/N in the thermo- 
dynamic limit N — > oo. To render this limit finite, the 
normalization 

/ N/2 \ -1 

Af N =\2j2<3)J (5) 

has been introduced in @. This normalization is a gen- 
eralization of the so-called Kac prescription commonly 
used when studying long-range interacting systems. For 



large system sizes A, its asymptotic behavior is of the 
form 

{(1 - a)2 1 - a N a ^ 1 for < a < 1, 
1/lnJV fora = l, (6) 

1/C(a) fora>l, 

where ( denotes the Riemann zeta function. Equation (|6[) 
implies that A/at vanishes in the thermodynamic limit for 
exponents a 1. 

In most studies of the (generalized) Ising model ©, 
observables related to the z components erf are consid- 
ered, like the magnetization per spin in the z direction, 
J2iLi a i/N- Since observables of this type commute with 
Hn, their expectation values are conserved quantities 
and equilibration cannot be observed. Here we want to 
consider the time-dependence of the expectation value of 
observables which are linear combinations of the x com- 
ponents erf of the spin operators, 

N 

A(ai,...,a N ) = ^2ai<jf (7) 

i=l 

with real coefficients a,. None of these observables com- 
mutes with the Hamiltonian ©. The set containing all 
A is, as pointed out in [14] . maximal Abelian, meaning 
that the information obtained by measuring simultane- 
ously all A cannot be improved by any compatible mea- 
surement performed at the same instant. Our aim is to 
study the time evolution of observables A, starting from 
an initial state (i.e., density operator) p(0) which is di- 
agonal in the erf tensor product eigenbasis of T-L. 

Related models and earlier results. — The above model 
is inspired by work of Emch [jjj and Radin [IB}. The 
crucial difference of our model, compared to the work 
of Emch and Radin, is the presence of strong long-range 
interactions and of the A-dependent normalization factor 
A/jv in the Hamiltonian ©, and we will see that these 
ingredients lead to dramatically different behavior. 

The model introduced in [14j, including the choices of 
observables and initial states, is inspired by induction de- 
cay experiments probing the pulsed magnetic resonance 
of nuclei in a CaF2 crystal [17J. In these experiments, 
the decay of the x component of the total magnetization 
A(l, . . . , 1) = YaLi a f was measured, and it was found 
to be superimposed by oscillations (see Fig. 1 of [171). 
Despite its simplifying assumptions, the model captures 
well both the decay and the beating observed experimen- 
tally. Moreover, the Emch-Radin model and its random 
generalizations [16( have proved useful as paradigmatic 
models for which the approach to equilibrium in quan- 
tum systems can be studied analytically. 

Time evolution. — We study the time evolution, gener- 
ated by H n ? of the expectation value of an observable A 
of the form ([7]) with respect to the initial state p(0), 

(A)(t) = Tr [e-' lHNt Ae iHNt p(0)] (8) 
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(in units where H = 1). Performing the trace in the erf 
eigenbasis of H, the diagonal form of p(0) implies that 
only the diagonal elements of the operator e _lifjvt ^4e lffjv * 
are required, and it is this crucial ingredient which allows 
us to obtain, similar to the calculation in (l4| . the exact 
result 

JV 

(A)(t) = (A)(0) cos(2W) Y[ cos 2 [2Af N e(j)t]. (9) 



In comparison with the original Emch-Radin model, the 
important difference here is the explicit N dependence 
of the argument of the cosine through the normalization 
jVjv. Regarding the approach to equilibrium, the Larmor 
precession cos(2fti) is not relevant and we set h = in the 
following. The behavior of ([5]) is plotted for exponents 
a = 2 and a = 1/2 and various system sizes A in Fig. 
[TJ and in all cases the expectation value of A appears 
to be decaying in time to the microcanonical ensemble 
average (A) mic — 0. However, this decay is only appar- 
ent, as we can read off from © that (A)(t) is an almost 
periodic function in time for all finite N, and Poincare re- 
currences will therefore occur on much longer time scales 
than shown. To possibly observe true equilibration, we 
have to invoke, as often in statistical physics, the idealiz- 
ing concept of the thermodynamic limit. In this limit, re- 
currence times may diverge (typically exponentially) and 
an approach to equilibrium may take place. An analysis 
of the infinite system dynamics is most rigorously done 
in a ^-algebraic language 1^, 16|, but it essentially boils 
down to discussing the large- A limit of the product in 
©. Similar to Lemma 4 of 
bound 



(A)(t)/(A)(0) 

1.0 



15| . we obtain the upper 



lim \{A)(t)\ ^ |{A)(0)| expl-GA/^ 2 } (10) 

where c is a positive constant. For a > 1, the large- N 
limit of the normalization Af x in ((6]) is strictly positive, 
proving a stretched exponential approach of (A) to its 
equilibrium value. For the case ^ a < 1 of strong long- 
range interactions we are particularly interested in, A/^ 
is zero. The bound in ([TU]) . therefore, gives the positive 
constant | (A) (0) | and fails to provide any indication on 
whether equilibrium is approached or not. Up to this 
point the analysis has been very much along the lines of 
earlier work reported in the literature. The main novel 
result of the present Letter is to complement the upper 
bound (|10p by a lower bound on (A) proving and char- 
acterizing the divergent equilibration time for ^5 a < 1 
in the large- TV limit. 

Proposition 1. Consider the Ising-type Hamiltonian 
Hn defined in ([2]) with power law decaying interactions 
e{j) = j~ a with ^ a < 1. Consider further an ob- 
servable A of the type ([7|) and an initial state p(0) being 
diagonal in the o~f tensor product eigenbasis of the un- 
derlying Hilbert space H. Then for the expectation value 
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FIG. 1. Time evolution of the expectation value ([9]) of an 
observable A for magnetic field h = and various system 
sizes N. Top panel: For a = 2, an apparent decay is observed, 
superimposed by oscillations. The N dependence of the time 
evolution is so weak that the curves for the various system 
sizes cannot be discerned in the plot. The behavior for other 
a > 1 is qualitatively similar. Bottom panel: For a — 1/2, 
the expectation value again appears to be decaying, but on a 
time scale that depends strongly on the system size iV (note 
the logarithmic scale). Similar behavior is observed for other 
values of a between zero and one. 



(A)(t) of A with respect to p{t), the following holds true: 
For any fixed time t and some small 8 > 0, there is a 
finite Nq(t) such that 

\(A)(t) - (A)(0)\ < 5 W < t, N > N (t). (11) 

Interpreting this result in terms of an experiment, we 
can think of an experimental resolution 5 for the mea- 
surement of A, and some duration r of the experiment. 
Then the above proposition states that, within the ex- 
perimental resolution and for a large enough system, no 
deviation of (A) (t) from its initial value can be observed 
for times t ^ r. In the above sense, (A)(t) converges in 
the thermodynamic limit to the constant (A)(0) which, 
in general, is different from the microcanonical ensemble 
average (A) mic = 0. 

Proposition [1] can be proved by constructing, for large 
enough A, a lower bound on the product in ^ by means 
of an integral approximation. The result can be ex- 
pressed in terms of hypergeometric functions and, by 
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virtue of the asymptotic properties of these functions, the 
bound can be pushed arbitrarily close to 1 by increasing 
the system size N. 

Comparing the scaling of the equilibration times tq 
with system size N, a striking difference between classi- 
cal and quantum mechanical quasistationary behavior is 
revealed: From the data plotted in Fig. Q] (bottom) and 
consistent with the bound (fT0| . one finds To oc N 1 with 
7 0.5, whereas in the classical case 7 ~ 1.7 has been 
extracted from numerical data and 7 = 1 from a Vlasov 
approach [18J. 

Generalizations. — To keep the presentation simple, 
the above exposition had been restricted to a one- 
dimensional lattice and power law interactions e(J) — 
j~ a , but several generalizations are straightforward. In 
fact, only the large-j asymptotic behavior of e(j) is rele- 
vant for the proof of Proposition [TJ and the result holds 
for any \e(j) \ decaying asymptotically as Cj~ a with some 
positive constant C and < a < 1. For lattice dimen- 
sions d, the results remain true for exponents < a < d, 
again with only marginal modifications of the proof. 

Instead of observables A of the type (|7J| we considered, 
linear combinations of y components of Pauli spin opera- 
tors would not have altered any of the conclusions. Linear 
combinations of z components, however, commute with 
the Hamiltonian (j2J and are therefore conserved quanti- 
ties which do not equilibrate. The results in [l5| indicate 
that even further generalizations should be feasible. 

Discussion and outlook. — We have analytically stud- 
ied the time evolution of long-range quantum spin models 
with Ising-type interactions © on lattices of arbitrary di- 
mension and for a certain class of initial states. For finite 
system sizes N, almost periodic behavior is observed with 
recurrence times that increase exponentially with N. In 
a rigorous sense, equilibration can therefore occur only in 
the thermodynamic limit of infinite lattice size. Whether 
this indeed happens depends on the asymptotic behavior 
of the interaction strength e(j) for large distances j on 



previously observed and extensively studied in classical 
systems, into the realm of quantum mechanics, while at 
the same time providing the first rigorous, analytic proof 
of quasistationary behavior from the microscopic time- 
evolution equations. 

It is tempting to speculate that quasistationary be- 
havior, i.e., equilibration times that diverge for large sys- 
tem sizes N, might show up under much more general 
conditions, or even generically for long-range interacting 
quantum systems. An extension of our results to general 
anisotropic Heisenberg models might be a promising fu- 
ture project and, although an exact solution for the time 
evolution of (A) seems to be out of reach, the derivation 
of suitable bounds on (A)(t) might be feasible. Finally, 
the simultaneous presence of short- and long-range inter- 
actions might lead to equilibration taking place on two 
different time scales, possibly leading to more complex 
behavior. The peculiarities of diverging relaxation times 
must be expected to be reflected also in many applied 
aspects of quantum spin systems, including the study of 
quantum quenches [2|, decoherence, and many others. 
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